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We compute the Gaussian-fluctuation corrections to the saddle-point Schwinger-boson results using 
collective coordinate methods. Concrete application to investigate the frustrated Ji — J2 antiferro- 
magnet on the square lattice shows that, unlike the saddle-point predictions, there is a quantum 
nonmagnetic phase for 0.53 < J2/J1 J$ 0.64. This result is obtained by considering the corrections to 
the spin stiffness on large lattices and extrapolating to the thermodynamic limit, which avoids the 
infinite-lattice infrared divergencies associated to Bose condensation. The very good agreement of 
our results with exact numerical values on finite clusters lends support to the calculational scheme 
employed. 



In the last years there has been a lot of interest in 
the properties of quantum magnetic systems, partic- 
ularly frustrated quantum antiferromagnets. Although 
this interest was initially related to the possible connec- 
tions between magnetism and superconductivity in the 
ceramic compounds, the current activity in the area is 
now well beyond this original motivation. 

Among the analytical methods used to study quan- 
tum spin systems, the Schwinger-boson approach |2j] is 
one of the most elegant and successful techniques. Con- 
trary to standard spin-wave theory, it does not rely on 
having a magnetized ground-state, which leads to nice 
rotational properties of the results and to the possibility 
of describing ordered and disordered phases in an unified 
treatment. However, this theory has the drawback of be- 
ing defined in a constrained bosonic space, with unphys- 
ical configurations being allowed when this constraint is 
treated as a soft (average) restriction. This drawback can 
be in principle corrected by including local fluctuations 
of the boson chemical potential. ||] 

Despite the widespread use in the literature of the 
Schwinger-boson representation of quantum spin opera- 
tors, we are not aware of a complete calculation of Gaus- 
sian corrections to saddle-point results. In particular, 
for frustrated quantum antiferromagnets such calcula- 
tions have been sketched several times, but never 
fully undertaken. In this work we fill up this gap by 
presenting the general calculation of Gaussian fluctua- 
tions in the Schwinger-boson approach. Since the theory 
presents a local U(l) symmetry, we use collective coordi- 
nate methods — as developed in the context of relativistic 
lattice gauge theories Q — to handle the infinitely-many 
zero modes associated to the local symmetry breaking in 
the saddle-point expansion. As a concrete application, 
we study the existence and location of the nonmagnetic 
quantum phase predicted to occur as a consequence of 
quantum fluctuations and frustration in the J1—J2 model. 

. . 

We will consider a general Heisenberg Hamiltonian, 



H = Yluj) JijSi-Sj, where (ij) are links on a lattice. 
We write spin operators in terms of Schwinger bosons: 
H Si = ^a\.a.ai, where at = (a].pa|^) is a bosonic spinor, 
o is the vector of Pauli matrices, and there is a boson- 
number restriction a^aia = 2S on each site. With 
this faithful representation of the spin algebra, the rota- 
tional invariant spin-spin interaction can be written as 
Si.Sj-.B^Bij : -AlAij. We defined the SU(2) invari- 
ants Aij = I J2a O a i<J a 3& and B ij = 5 J2a a \a a 3<y = 

—a, a = ±), and the notation :0: indicates the normal 
order of operator O. 

Using boson coherent states to formulate the par- 
tition function, we formally integrate the Schwinger 
bosons by decoupling the quartic terms in H using two 
complex Hubbard-Stratonovich fields W[j (r) (r = 1,2 
correspond, respectively, to the fields that couple to 
the ferromagnetic B and antiferromagnetic A channels). 
This leads to Z = J[VWVW] [V\] e - A , where the mea- 
sure [VWVW]=Y\ (lj) ^M 3 /^)dW 1 i0 {r)dWl 3 (r) and 

[PA] = Yli T d\i(T) /2tt. Here W is the complex conju- 
gate of W and the effective action 

\m,r i J 

+T h tW,W,\). (I) 

The integration on A comes from the integral represen- 
tation of the ^-functions which force the boson number 
restriction on each site. In (|l]) T\,{W ,W,\) is the free 
energy of a boson gas coupled to W, W. In a compact 
notation where site (i), Nambu (s) and Matsubara (r) 
indices are summed over in the trace (Tr) operation, we 
have F h (W, W, A) = Tr \nM(W, W, A). The dynamical 
matrix Adfj (r) is given by 
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where we have suppressed the r dependences to simplify 
notation. It is important to stress that this formal result 
for J-b must be understood as the limit of a discrete-r 
mesh. ]8| In the functional formulation such a procedure 
picks up the zero-point contributions to the energy com- 
ing from the Bogoliubov transformation in the canonical 
approach. 

The action A is invariant under both local gauge trans- 
formations, WT- -> WT.e^ i±( ^ (- for r = 1, + for r = 2), 
Xi — > \i—d r 8i, and global SU(2) rotations of the original 
bosonic spinors. Consequently, it would be convenient to 
include external sources fj, r\ which explicitly break the 
global SU(2) invariance: T\> — > T\> — r\Mr x r\. They are 
useful to study the spontaneous breaking of this sym- 
metry at T — as a consequence of Bose condensation, 
which signals the onset of magnetic long-range order in 
the ground state. || However, in order to avoid the in- 
frared divergencies associated to this symmetry breaking, 
in what follows we approach the thermodynamic limit by 
extrapolation of results on large but finite lattices. 

Up to this point all the manipulations are exact. To 
proceed further we resort, as usual, to a saddle-point ex- 
pansion. Using the collective notation for the fields <p< = 
(W, W, A), the (static) saddle-point values W lj ,W[ j 
and A^o are obtained from the extremal equations 
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where $0 = ( JW , JW , 2Si) and Go = M 1 is the 
bosonic Green function at saddle-point order. || After 
expanding A to second order in the field fluctuations 
around these saddle-point values, we end up with 
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In this equation Aq = PJ2{ij) J ij (l^ij'.ol 
the effective action evaluated at the saddle point, and 
A0 are the field fluctuations. The matrix A^ is given 
by 
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where J is a diagonal matrix containing the exchange 
couplings Jij along the diagonal, except for the entries 
corresponding to A — A derivatives that are zero. 

We stress that the saddle-point expansion of the ef- 
fective action breaks the gauge invariance of the the- 
ory. Although the spontaneous breaking of a local sym- 
metry is strictly forbidden by Elitzur's theorem, a 
possible justification of this procedure has been given 
in the context of relativistic lattice gauge theories: Jll]] 
The existence of saddle-point solutions connected by the 



continuous U(l) gauge group makes the A^ 1 matrix to 
have infinitely-many zero modes, which are the Gold- 
stone bosons associated to the spurious local symme- 
try breaking. In particular, for translational-invariant 
saddle-point values W[ — W[_j , A^ = Ao, trans- 
forming to momentum-frequency space there is a zero 
mode ipg(k,oj) — S(p(d)/S9] iU j\sp in every k— uj subspace. 
Here 0(9) is the vector of gauge-transformed fields (No- 
tice that (f>o' is a right eigenvector of the nonhermitian 
matrix „4( 2 )). To avoid the infinities associated to these 
modes without restoring forces — which correspond to lo- 
cal symmetry transformations — we introduce collective 
coordinates along the gauge orbit. Exact integration of 
these coordinates eliminates the zero modes and restore 
the gauge symmetry (in the sense that noninvariant oper- 
ators average to zero). JlT| This program can be carried 
out by enforcing in ||) the so-called background gauge 
condition or "natural" gauge, || 0(6)^ .0$ = 0. This 
condition can be introduced into the functional measure 
by means of the Faddeev-Popov trick, and restricts the 
integration to fields fluctuations which are orthogonal to 
the collective coordinates. At T = 0, after carrying out 
the remaining integrations in (^|) on the genuine fluctu- 
ations, we obtain the one- loop correction to the ground- 
state energy per site, 
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Here the Fadeev-Popov determinant App(k,o;) = 
\0q (k, Lu).0Q(\t, oj)\ (0q(\s.,uj) is the left zero mode of 

_4.( 2 ) in the k— u> subspace), and A^ is the projection of 
„4^ 2 ) in the subspace orthogonal to the collective coordi- 
nates. 

We have particularized these results for a matrix 
which couples a given site to its first (Ji) and second 
( J2) neighbors on a square lattice. This model with frus- 
trating Ji , J2 > has been much studied in the last 
years in connection with the physics of the lightly doped 
Cu02 planes (see |l2j and references therein). The main 
problems under discussion are the location and physical 
properties of a possible nonclassical phase for intermedi- 
ate frustration a = J2/J1 ~ 1/2. In the ordered phases 
of this model (with magnetic wavevectors Q = (jr, it) for 
small a and Q = (w, 0) for a <~ 1) the Bose condensate 
breaks the global SU (2) symmetry and its density gives 
the local magnetization. || The associated physical Gold- 
stone modes at k = 0, Q lead to serious infrared divergen- 
cies of intermediate quantities, which have to be cured 
by standard renormalization prescriptions. In order to 
avoid these problems we have computed physical quan- 
tities (which are free of these divergencies) on large by 
finite lattices, and finally extrapolated these values to the 
thermodynamic limit. However, this procedure requires 
a careful numerical treatment, since finite results are ob- 
tained by cancellations among products of the very large 
entries of the matrix A^ . 
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The equations (||)-(||) lead to a ground-state energy 
Eqs =Eq + Ei as shown in Fig. 1. This figure contains 
the result for the infinite lattice and also for finite lattices 
of 16 and 36 sites, which allow a comparison with exact 
results obtained by numerical diagonalization. |l3| ] The 
addition of the Gaussian correction E\ (eq. (Q)) slightly 
improves the already very good saddle-point value Eq- At 
saddle-point order the theory predicts a first order tran- 
sition between the two magnetic ground states at some 
intermediate frustration (a ~ 0.6), with no intervening 
disordered phase. This is in contradiction with numer- 
ical studies and series expansion results, [H| which 
found a disordered purely-quantum phase somewhere in 
the range 0.3 < a < 0.6. However, it has been recently 
shown that these finite-lattice studies probably did not 
reach the scaling region where the required extrapola- 
tion can be trusted. |16| In our case, the existence or not 
of magnetic long-range order was investigated by con- 
sidering the spin-stiffness tensor p s — d 2 i?(Q)/<9Q<9Q, 
where £gs(Q) I s the ground-state energy with twisted 
boundary conditions. Qjl^l In particular, the correc- 
tions to the saddle-point values obtained in |l(| lead to 
the results shown in Fig. 2. It gives the stiffness along 
one of the lattice directions in the Q = (tt, tt) phase at 
saddle-point and one-loop order, on lattices of 16 and 20 
sites (according to numerical studies H] the tilted 20- 
site lattice has a stiffness closer to the thermodynamic 
limit than regular clusters). The comparison with exact 
results [[l3| shows a good agreement, especially for the 
20-site lattice. The extrapolation to the thermodynamic 
limit leaves a window 0.53 < a < 0.61 where the stiff- 
ness vanishes and the magnetic order is melted by the 
combined action of quantum fluctuations and frustration 
(see Fig. 3). This result is in fairly good agreement with 
the second-order modified spin wave calculation of [ pS] , 
which predicts a quantum disordered phase in the range 
(0.52, 0.57). Notice, however, that our theory -valid also 
in the disordered region- predicts that for 0.61 < a < 0.64 
the short-range antiferromagnetic order is energetically 
more favourable than the weak collinear order (see dotted 
line in Fig. 1). Finally, we have computed the spin- wave 
velocity c s in the Neel phase using the finite-size formula 
E GS (N) ~ E GS - a/N 3 / 2 and the relation a = 1.4372c s . 
Jl8[ For J 2 = we obtained c s ~ 1.52 Ji at saddle-point or- 
der, and c s c±1.37Ji including the Gaussian fluctuations. 
These values should be compared to the result of exact 
diagonalization, |ll| c s ~ 1.28 — 1.44 Ji, which depends on 
the cluster sizes considered to extrapolate. In Fig. 3 we 
plot the scaling coefficient a as a function of frustration 
(in the collinear phase this coefficient should be propor- 
tional to some anisotropy-averaged velocity) . In the Neel 
phase the spin- wave velocity vanishes at a ~ 0.54, slightly 
above the point where the stiffness goes to zero. In the 
collinear phase the averaged velocity never vanishes, and 
actually it blows up at a ~ 0.58. Notably, this is the value 
where the stronger stiffness in the antiferromagnetic di- 
rection of the collinear order vanishes. More details and a 
further discussion of these results will be given elsewhere. 



Before closing, it is worth mentioning a few side ques- 
tions we have considered while performing these calcu- 
lations. First, since the Hamiltonian commutes with 
the local boson- number restrictions, we considered im- 
posing them only through r-independent Lagrange mul- 
tipliers Xi. This reduces the local gauge symmetry 
to static transformations and, consequently, A^ con- 
tains zero modes only in the lu = subspace. For 
the calculation of the ground-state energy — which in- 
volves an integration over u> — it is then allowed to 
simply forget about these modes and compute E\ = 

47/-!L dw Ek lndct A ( runc( k >^)> wh ere ^runc is ec l ual 

to A^ with the column and row corresponding to the 
(static) A fluctuations deleted. In this case there is no 
need for App and (|J) is equivalent to the RPA result for 
the bosonic theory with Ao taken as a chemical potential. 
We checked that this procedure gives the same result that 
(||) , although the numerical evaluation requires extra care 
because of the divergencies near w = 0. Second, in most 
works in the literature on Schwinger bosons the identity 
: B\-Bij : +A\-A^ = S 2 (which holds because of the con- 
straint) is used to simplify the Hamiltonian, leaving only 
the ferromagnetic B or the antiferromagnetic A chan- 
nel in the formulation. In a previous publication |l9| we 
warned that at saddle-point order this produces large er- 
rors in the ground-state energy, since the above identity 
is largely violated when the constraint is imposed only on 
average. We proved that this remains true even after the 
inclusion of fluctuations, with the contributions of both 
channels being important. Third, in (|j) we can write 

In det A ( l } = Tt In A^ = - £^=1 ±Tr{I - Af) 11 . If only 
the first term (n = l) is added to the saddle-point energy 
Eq, the result is equivalent to the energy obtained in a 
full Hartree-Fock-Bogoliubov decoupling of the quartic 
boson interactions in H . Moreover, it can be proved that 
Ehfb = Eq + E n= i = §i?o- That is, a fully self-consistent 
treatment takes advantage of the unphysical enlargement 
of Fock space to lower the energy too much. The inclu- 
sion of terms with n > 2 corrects -Ehfb in nearly 50%. 

In conclusion, we have computed the Gaussian- 
fluctuation corrections to the Schwinger-boson saddle- 
point results using collective coordinate methods. As 
a concrete application, we investigated the ground-state 
structure of the J\ — Ji model. By considering the spin 
stiffness of this model we showed that, contrary to the 
saddle-point predictions, there is a quantum nonmag- 
netic phase that intervenes between the Neel (tt, tt) and 
collinear (n, 0) orders, as suggested by numerical meth- 
ods. Its stability region is, however, smaller than pre- 
dicted by these methods. Moreover, the comparison with 
exact results on finite lattices lends support to our cal- 
culations, which, in addition, have the advantage of a 
well-defined thermodynamic limit. Finally, the theory 
developed here can be easily extended to spiral phases, 
]l9t which would allow to investigate, for instance, the 
ground-state order in triangular |20| and kagome lattices. 
pl| Work in this direction is in progress. 
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FIG. 1: Ground-state energy per site Eqs as a func- 
tion of J2/J1 for lattices of (from top to bottom) N = 00, 
36, and 16 sites. The full and dashed lines give the 
fluctuation-corrected and saddle-point results, respec- 
tively. Points are exact results from [fl3|| . The dotted 
line in the range 0.53 < J2/J1 ^ 0.64 indicates the re- 
gion without long-range magnetic order. Notice that for 
clarity the curves corresponding to N — 00 and 36 have 
been shifted upwards in 0.1 and 0.05 respectively. 



FIG. 2: Spin stiffness p s in the Neel phase for a) a 16- 
site lattice, and b) a tilted 20-site lattice. Symbols and 
linotypes are the same as in Fig. 1. 

FIG. 3: Scaling coefficient a proportional to the spin- 
wave velocity (top panel) and spin stiffness p s (bottom 
panel) extrapolated to the thermodynamic limit. Line- 
types are the same as in Fig. 1. 
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